REMARKS 

Claims 1-41 are pending in this application and stand rejected. No claims have 
been amended in this response and no claims have been added. The specification has 
been amended to add paragraph numbers. No new matter has been added. Applicant 
respectfully requests reconsideration of this application and its pending claims in light of 
the following remarks. 

I. Drawings 

Applicant respectfully submits formal drawings (Figures 1-43) herewith in 
compliance with 37 CFR 1.84(f). 

II. Request for Information 

Applicant respectfully submits of the following publications as cited in the instant 

application for the Examiner's review. 

" Fundamentals of Geostatics in Five Lessons ", Journel, A.G., 
Short Course in Geology, Vol. 8, 44 pp. 15-16, AGU, Washington, D.C., 
1989. 

" GSLIB Geostatistical Software Library and User's Guide second 
edition ", Deutsch, Clayton V. and Journel, A., Oxford University Press, 
pp. 14-15, New York, Oxford, 1998. 

III. The Specification 

A substitute specification is submitted herewith, amended to include paragraph 
numbering. No amendments have been made. No new matter has been added. 

IV. The § 103 Rejections 

In the Office Action dated April 8, 2002, the Examiner rejected claims 1-41 under 
35 USC 103(a) as being unpatentable. 



94.0027 



14 



A. Claim 1 



Claim 1 is rejected under 35 U.S.C. 103(a) as being unpatentable over Jones in 

view of Matteucci and Tucker. Applicant respectfully traverses this rejection. 

Claim 1 recites: 

A method of generating a map illustrating a set of 
characteristics of a cross section through an earth formation 
representing a time slice or a horizon through said 
formation in response to a plurality of scattered data 
observations on said cross section representing a plurality 
of parameters located at a plurality of locations on said 
cross section, comprising the steps of: 

(a) gridding said cross section thereby generating a gridded 
cross section which includes a grid having a plurality of 
intersections and said plurality of scattered data 
observations distributed among the intersections of said 
grid on said cross section; 

(b) obtaining a unique cumulative distribution function 
associated with each intersection of the grid of the gridded 
cross section thereby producing a plurality of cumulative 
distribution functions associated, respectively, with the 
plurality of intersections of said grid; 

(c) choosing a value from each of the cumulative 
distribution function at each of the intersections of the 
gridded cross section thereby producing a plurality of 
values associated, respectively, with the plurality of 
intersections, and 

(d) assigning each value to its associated intersection of the 
gridded cross section and assigning a unique color to said 
each value thereby generating a map illustrating said set of 
characteristics of said cross section through said earth 
formation. 



The Examiner states that Jones (U.S. Patent No. 5,838,634) discloses "a gridded 
cross section" with its disclosure of a three dimensional array of cells and mention of 
"structural surfaces or horizons in the form of 2-D computer grids or meshes." Applicant 
agrees that Jones discloses 3-D arrays of cells, which are known in the art and the 3D 
array inherently includes cross sections of horizons. As the Examiner points out, Jones 
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mentions modeling horizons or structural surfaces with meshes or grids. Applicant notes 
that Jones does not mention modeling time slices with grids and Jones is not directed to 
any further analysis with respect to gridded cross sections. Instead, Jones divides the 
formation up into blocks (not cross sections) and the modeling described by Jones is 
performed on a block basis, see Jones column 12, lines 46 through column 13, line 29, 
not through the use of cross sections. Although grids by definition have intersections, 
Jones thus does not disclose "said plurality of scattered data observations distributed 
among the intersections of said grid on said cross section... " or the other limitations of 
claim 1 , as acknowledged by the Examiner. 

Matteucci (U.S. Patent No. 5,884,229) "is a statistical method for quantitatively 
measuring the lateral continuity of the seismic reflection character at in specified location 
in a subsurface target formation." (Mateucci, column 3, lines 35-38.) In Mateucci, the 
analysis is focused on a comparison of seismic traces, not 3-D arrays or gridded cross 
sections. A seismic data trace is selected as a "reference trace." Note that each seismic 
trace is a response to sound waves bouncing off underground formations and can be seen 
as a vertical , wavy line 20, passing through the formation in the z direction in Mateucci, 
FIG. 1. In Mateucci, "all [seismic] traces falling within a specified distance from the 
reference trace are extracted from the survey, and one or more statistics which measure 
the similarity between the reference trace and each of the other traces are calculated. . . . 
Finally, the calculated statistics are used to make a determination of whether or not the 
reflection character within the target interval of each of the extracted traces is the same as 
or different from the reference trace." (Mateucci, column 5, lines 29-44). While 
Mateucci uses the cumulative distribution function statistical tool, it is within the context 
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of this comparison of seismic traces and each cumulative distribution function is taken 
for a particular seismic trace (see Mateucci column 7, lines 2-15). This is a wholly 
different context than that of the present invention, where the cumulative distribution 
functions are taken at each intersection of the grid on the cross section . Mateucci goes 
on to take the Kolmogorov-Smirnov statistic, which "measures the similarity of the 
amplitude distribution of the two seismic traces." (Column 6 line 58- 60 and continuing 
through column 7, line 15). Again, Mateucci is using the cumulative distribution 
functions to measure the similarity of two vertical seismic traces, not investigating 
similarity of points along a grid on a cross section through the formation. 

If one were to combine Jones and Mateucci, the result would be a 3-D cube array 
or block, having horizons gridded but with no further analysis of the gridded horizons. 
The array would have a plurality of seismic traces passing vertically through it and the 
invention would involve an analysis of the similarity of seismic traces using cumulative 
distribution functions and Kolmogorov-Smirnov statistics. This does not meet the 
limitations as recited by claim 1, which include "a plurality of scattered data observations 
distributed among the intersections of said grid on said cross section; obtaining a unique 
cumulative distribution function associated with each intersection of the grid . . . choosing 
a value from each of the cumulative distribution function at each of the intersections of 
the gridded cross section . . . and assigning each value to its associated intersection of the 
gridded cross section and assigning a unique color to said each value ..." (Emphasis 
added.) 

Nor does Tucker ( The Computer Science and Engineering Handbook , edited by 
Allen B. Tucker, Jr., 1997) supply the deficiencies of Jones and Mateucci. While 
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Tucker does disclose the use of color in "a particular slice of 3-D dataset," it does not 
disclose the use of color as recited in claim 1, i.e., assigned to a value taken from the 
cumulative distribution function of an intersection on the grid of the cross section. 

Accordingly, the combination of Jones, Mateucci and Tucker does not disclose 
"a plurality of scattered data observations distributed among the intersections of said grid 
on said cross section; obtaining a unique cumulative distribution function associated with 
each intersection of the grid; ... choosing a value from each of the cumulative 
distribution function at each of the intersections of the gridded cross section; ... and 
assigning each value to its associated intersection of the gridded cross section and 
assigning a unique color to said each value" as recited by claim 1 (emphasis added). Thus 
claim 1 as amended is felt to distinguish patentably over the combination of the Jones, 
Mateucci and Tucker references. 
B. Claim 2 

Claim 2 is rejected under 35 U.S.C. 103(a) as being unpatentable over Jones in 
view of Matteucci and Tucker and Journel (Fundamentals of Geostatics in Vive 
Lessons," by Journel, vol. 8AGU, 1989). Applicant respectfully traverses this rejection. 

As the Examiner mentions, claim 2 depends from claim 1 with three 
additional limitations, which are not disclosed by Jones: (bl) Krigging the gridded cross 
section thereby generating a plurality of expected values and a plurality of corresponding 
standard deviations associated, respectively, with the plurality of intersections of the grid 
of the gridded cross section; (b2) producing a probability density function associated with 
each expected value and each corresponding standard deviation generated from step (bl) 
thereby producing a plurality of probability density functions corresponding, 
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respectively, to the plurality of intersections of the grid of the gridded cross section; and 
(b3) producing a cumulative distribution function associated with each probability 
density function produced from step (b2) thereby producing a plurality of cumulative 
distribution functions corresponding, respectively, to the plurality of probability density 
functions which correspond, respectively, to the plurality of intersections of the grid of 
the gridded cross section. 

Claim 2 is dependent on claim 1 and contains all of its limitations so it, too, is felt 
to be distinguishable from the cited references. With respect to its additional limitations, 
the Examiner states that Krigging and probability functions are described in the instant 
Specification, which is correct, but the use of such tools as recited by claim 2 in 
connection with being taken at and corresponding to intersections of the gridded cross 
section as recited by claim 2. Similarly, as previously pointed out, Matteucci discloses 
the use of cumulative distribution functions, but for use in comparing seismic traces, not 
taken at intersections of grids on gridded cross sections. Accordingly, claim 2 is felt to 
be patentably distinguishable from the cited references. 
C. Claim 3 

Claim 3 is rejected under 35 U.S.C. 103(a) as being unpatentable over Jones in 
view of Matteucci and Tucker and Journel and Hogg (Probability and Statistical 
Inference, by Hogg, Robert V and Tanis, Eliot A, 3d ed. 1988). Claim 3 is dependent on 
claim 1 and contains all of its limitations so it, too, is felt to be distinguishable from the 
cited references. The addition of Hogg does not affect this as the reference does not 
disclose the application of the statistical tool of choosing a probability from the 
cumulative distribution function in the circumstance where the cumulative distribution 
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function is taken at an intersection of a grid on a gridded cross section. Such an 
application is not disclosed by Hogg, which does not supply the deficiencies of the 
previously discussed references. Accordingly, Applicant respectfully traverses this 
rejection. 

D. Claim 4 

Claim 4 is rejected under 35 U.S.C. 103(a) as being unpatentable over Jones in 
view of Matteucci and Tucker and Journel and Hogg. Applicant respectfully traverses 
this rejection. Claim 4 depends from claim 2, contains all its limitations and is felt to be 
likewise patentable over the cited references. As described with respect to claim 3, the 
addition of Hogg does not supply the deficiencies of the previously discussed references. 
Therefore, Applicant respectfully asserts it has traversed this rejection. 

E. Claim 5 

Claim 5 is rejected under 35 U.S.C. 103(a) as being unpatentable over Jones in 
view of Matteucci and Tucker and Journel and Hogg. Claim 5 depends from claim 2, 
contains all its limitations and is felt to be likewise patentable over the cited references. 
As described with respect to claim 3, the addition of Hogg does not supply the 
deficiencies of the previously discussed references. Therefore, Applicant respectfully 
asserts it has traversed this rejection. 

F. Claim 6 

Claim 6 is rejected under 35 U.S.C. 103(a) as being unpatentable over Jones in 
view of Matteucci and Tucker and Journel and Hogg. Claim 6 depends from claim 2, 
contains all its limitations and is felt to be likewise patentable over the cited references. 
As described with respect to claim 3, the addition of Hogg does not supply the 
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deficiencies of the previously discussed references. Therefore, Applicant respectfully 
asserts it has traversed this rejection. 

G. Claim 7 

Claim 7 is rejected under 35 U.S.C. 103(a) as being unpatentable over Jones in 
view of Matteucci and Tucker and Journel and Hogg. Claim 7 depends from claim 2, 
contains all its limitations and is felt to be likewise patentable over the cited references. 
As described with respect to claim 3, the addition of Hogg does not supply the 
deficiencies of the previously discussed references. Therefore, Applicant respectfully 
asserts it has traversed this rejection. 

H. Claim 8 

Claim 8 is rejected under 35 U.S.C. 103(a) as being unpatentable over Jones in 
view of Matteucci and Tucker and Journel and Webber (U.S. Patent No. 6,081,577). 
Claim 8 depends from claim 2 and contains all of its limitations, so it is felt to be likewise 
patentable over the cited references. Webber does not supply the deficiencies of the 
previously discussed references. While it does, as the Examiner points out, disclose the 
use of affine correction, it is in the context of the transformation of three dimensional 
matrix volumes and des not disclose "applying an affine correction to each of the values 
chosen from each of the cumulative distribution functions associated with each of the 
intersections of the gridded cross section" as recited by claim 8. Therefore, Applicant 
respectfully asserts that this rejection is traversed. 

I. Claim 9 

Claim 9 is rejected under 35 U.S.C. 103(a) as being unpatentable over Jones in 
view of Matteucci and Tucker and Journel and Webber. Claim 9 depends from claim 
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8 and contains all of its limitations, so it is felt to be likewise patentable over the cited 
references. Therefore, Applicant respectfully asserts that this rejection is traversed. 
J. Claims 10- 14 

Claims 10-14 are rejected under 35 U.S.C. 103(a) as being unpatentable over 
Jones in view of Matteuci and Tucker and Journel and Hogg. Claims 10-14 depend 
from claim 9, contain all of its limitations, and so are felt to be likewise patentable over 
the cited references. Therefore, Applicant respectfully asserts that this rejection is 
traversed. 
K. Claim 15 

Claim 15 is rejected under 35 U.S.C. 103(a) as being unpatentable over Jones in 
view of Matteucci and Tucker, for the reasons described with respect to claim 1 . Claim 
15 is an independent "Program storage device ... cross section claim" with limitations 
corresponding to those of claim 1. Therefore, claim 15 is felt to be patentably 
distinguishable from the cited references for the same reasons discussed above with 
respect to claim 1, i.e. that the cited references do not disclose or suggest "obtaining a 
unique cumulative distribution function associated with each intersection of the grid of 
the gridded cross section thereby producing a plurality of cumulative distribution 
functions associated, respectively, with the plurality of intersections of said grid" or 
"choosing a value from each of the cumulative distribution function at each of the 
intersections of the gridded cross section ...;" or "assigning each value to its associated 
intersection of the gridded cross section and assigning a unique color to said each 
value. . . ." Therefore, Applicant respectfully asserts that this rejection is traversed. 
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L. Claims 16-24 

Claims 16-24 are rejected under 35 U.S.C. 103(a) as being unpatentable over 
Jones in view of Matteucci and Tucker, and, variously, in combination with Journel, 
Hogg and Webber for the reasons described with respect to claims 2-14. Claims 16-24 
are dependent, directly or indirectly, on claim 15 and contain all of its limitations. For 
the reasons described above with respect to claim 15 and claims 2-14, it is believed that 
claims 16-24 are patentably distinguishable from the cited references and that this 
rejection is traversed. 
M. Claim 25 

Claim 25 is rejected under 35 U.S.C. 103(a) as being unpatentable over Jones 
in view of Matteucci and Tucker. Claim 25 is an independent apparatus claim with 
limitations corresponding to those of claim 1. Claim 25 is felt to be patentably 
distinguishable from the cited references for the same reasons discussed above with 
respect to claim 1, i.e. that the cited references do not disclose or suggest " a second 
apparatus responsive to said first gridded cross section adapted for Kriging said first 
gridded cross section thereby generating a second gridded cross section having a 
plurality of intersections wherein each intersection of said second gridded cross 
section includes an expected value of a parameter and a standard deviation ; third 
apparatus responsive to said second gridded cross section for generating a plurality of 
cumulative distribution functions associated, respectively, with said plurality of 
intersections of said second gridded cross section ; and fourth apparatus adapted for 
selecting a plurality of values, respectively, from said plurality of cumulative 
distribution functions and for assigning said plurality of values and a plurality of 
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unique colors to the respective plurality of intersections of said second gridded cross 
section thereby generating said map. 55 Therefore, Applicant respectfully asserts that 
this rejection is traversed. 
N. Claims 26-29 

Claims 26-29 are rejected under 35 U.S.C. 103(a) as being unpatentable over 
Jones in view of Matteucci and Tucker and, variously, in combination with Journel, 
Hogg and Webber for the reasons described with respect to claims 2-14. Claims 26-29 
are dependent, directly or indirectly, on claim 25 and contain all of its limitations. For 
the reasons described above with respect to claim 25 and claims 2-14, it is believed that 
claims 26-29 are patentably distinguishable from the cited references and that this 
rejection is traversed. 
O. Claim 30 

Claim 30 is rejected under 35 U.S.C. 103(a) as being unpatentable over Jones in 
view of Matteucci and Tucker for the reasons described with respect to claim 1 . Claim 
30 recites: 

A method of generating a cube illustrating a set of 
characteristics of an earth formation disposed within a 
cubic volume of earth, said cube including a plurality of 
cross sections, each cross section including a plurality of 
scattered data samples, each cross section being gridded 
and including a plurality of intersections, comprising the 
steps of: 

(a) determining a plurality of cumulative distribution 
functions corresponding, respectively, to the plurality of 
intersections for each of said plurality of cross sections; 

(b) selecting a value from each of said cumulative 
distribution functions thereby selecting a plurality of values 
corresponding, respectively, to said plurality of cumulative 
distribution functions for each of said plurality of cross 
sections; 
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(c) assigning said plurality of values, respectively, to said 
plurality of intersections for each of said plurality of cross 
sections; and 

(d) assigning a plurality of unique colors, respectively, to 
said plurality of values assigned, respectively, to said 
plurality of intersections. 

As described in this response to claim 1, the cited references do not disclose 
"determining a plurality of cumulative distribution functions corresponding, 
respectively, to the plurality of intersections for each of said plurality of cross 
sections ," or "selecting a value from each of said cumulative distribution functions 
thereby selecting a plurality of values corresponding, respectively, to said plurality of 
cumulative distribution functions for each of said plurality of cross sections;" or 
"assigning said plurality of values, respectively, to said plurality of intersections for 
each of said plurality of cross sections," or "assigning a plurality of unique colors, 
respectively, to said plurality of values assigned, respectively, to said plurality of 
intersections." Therefore, Applicant respectfully asserts that this rejection is 
traversed. 
P. Claims 31-35 

Claims 31-35 are rejected under 35 U.S.C. 103(a) as being unpatentable over 
Jones in view of Matteucci and Tucker and, variously, in combination with Journel, 
Hogg and Webber for the reasons described with respect to claims 2-14. Claims 31-35 
are dependent, directly or indirectly, on claim 30 and contain all of its limitations. For 
the reasons described above with respect to claim 30 and claims 2-14, it is believed that 
claims 31-35 are patentably distinguishable from the cited references and that this 
rejection is traversed. 
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Q. Claim 36 

Claim 36 is rejected under 35 U.S.C. 103(a) as being unpatentable over Jones in 
view of Matteucci and Tucker. Claim 36 is directed to a "program storage device 
readable by a machine, tangibly embodying a program of instructions executable by the 
machine to perform method steps for generating a cube..." and stands rejected for the 
reasons described with respect to claim 1 . 

As described in this response to the rejections cited against claim 1, the cited 
references do not disclose "determining a plurality of cumulative distribution functions 
corresponding, respectively, to the plurality of intersections for each of said plurality of 
cross sections;" or "selecting a value from each of said cumulative distribution functions 
thereby selecting a plurality of values corresponding, respectively, to said plurality of 
cumulative distribution functions for each of said plurality of cross sections;" or 
"assigning said plurality of values, respectively, to said plurality of intersections for each 
of said plurality of cross sections;" or "assigning a plurality of unique colors, 
respectively, to said plurality of values assigned, respectively, to said plurality of 
intersections." Therefore, Applicant respectfully asserts that this rejection is traversed. 
R. Claims 37-41 

Claims 37-41 are rejected under 35 U.S.C. 103(a) as being unpatentable over 
Jones in view of Matteucci and Tucker and, variously, in combination with Journel, 
Hogg and Webber for the reasons described with respect to claims 2-14. Claims 37-41 
are dependent, directly or indirectly, on claim 36 and contain all of its limitations. For 
the reasons described above with respect to claim 36 and claims 2-14, it is believed that 
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claims 37-41 are patentably distinguishable from the cited references and that this 
rejection is traversed. 
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• 



CONCLUSION 



It is respectfully submitted that this application, as now amended, is in condition 
for allowance for the reasons stated above. Applicant respectfully requests 
reconsideration of this application and allowance of all its pending claims 

This amendment is intended to be a complete response to the Office Action dated 
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(30) 



with: 



2-(* < ,) a3 i 0 + J^ lsl 6 fl /(x # ~x o )= datum z(i„) 

for aD q = l,-",n 
(31) 

Then one may argue on the arbitrariness of choosing, in 
the SK case, covariance functions for interpolation func- 
tions. As was shown before, the choice of eovariance 
functions allows to minimize the resulting error variance, 
Var{Z 0 -£ 0 -}. 
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Ill: Linear Regression Under Constraint(s) and 
Ordinary Kriging (OK) 



In all the simple kriging (SK) developments of the 
previous Lesson II, the mean(s) were supposedly known, 
cf. relations (9) and (18). Consider, for example, the; 
SK estimator (18) of an unknown Z(z Q ) from n data 
related to the same attribute but at different locations 
SaiOt s l y * - *,n: - 



(1) 



The common (stationary) mean m is supposed known. 
If the process Z(x) was sampled at all locations x within 
the formation A considered, that mean m would be known 
but also there would not be any estimation problem left. 
Thus, in practice, either m is to be estimated prior to 
the SK algorithm, or an estimation algorithm should be 
designed which does not require prior estimate of that 
mean. 

Consider the linear estimator: 

ar=l 

The error mean is: 

E{Z{x 0 )-Z m (x 0 )} = m-ko-jr X a m = -A 0 +m fl - £ A a ) 

ami \ a J 



and should be set to zero, whatever the unknown value of 
m. This can be achieved only if: 

Y " / A o = 0 

Thus, an unbiased linear estimator of Z(xo) is written as: 

Z m (*o) » EAo£(z 0 ) , with: £^ = 1 (2) 

Note. Although for reason of convenience, we have used 
the. same notation X 0 for the SK weights of relation (l) 
and the weights of relation (2) they are not equal. 
" The error variance is written, as in Lesson II (11): 

a\ = Vnr {Z(* 0 )-Z-(xo)} « £ £ (3) 

cr=O0=Q 

with:!^; 1 ; 

and: C{x { 



x 0 )=Cov{Z(x a ),Z(z 0 )} 



We wish to minimize that error variance, still ensur- 
ing unbiasedness f the estimator, i. ensuring the con- 
straint: Xa = 1- This amounts to an optimization 
under linear constraint. The Lagrange formalism will 
be used whereby a function, of the (n+ 1) parameters 
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KuCt »V--,nand 2p bong lb Lagrange parameter, is 
defined as: 

-£(Ai,a-l,. .•<,!»; /i) = a % + 2 * Ao -lj 

The extreme of that function S is obtained by setting to 
zero its (n + l) partial derivatives: 

= ^0=1 A B — 1 ^ fl 

Th n first equations are rewritten: 

n 

C(x a - 1 0 ) - £ - **) - /* - 0 , i.e., 

0=1 

■ 

]£ - X fl ) + p « C(* 0 - 1 0 ) , a = 1, • • • ,1 



0=1 



Finally, the (n+1) unknowns, the V« and /*, are given by 
a system of (n + 1) linear equations, known as the "eon- 
strained normal" system of equations, belatedly renamed 
"ordinary krigmg" (OK) system: 

!£JLi -V?(*o -*/>) + 11 = C(* 0 -i 0 ) , ot = l,-..,n 

^ (4) 
The corresponding minimized error variance, ca l led 
"ordinary kriging" (OK) variance is written: 

0&jr = -3{[-?(*o)-^(* 0 )] 2 } = 



» C(0) - A B C(x a - s„) - /i > 0 



(5) 



Just like the SK system of Lesson II (12), the OK 
system (4) presents one and only one solution as soon as 
the data covariance matrix K = [C(x a - x 5 )] is positive 
definite, in practice as soon as: 

(i) the covariance function C(h) is licit, Le., is a positive 
definite function, see Lesson I condition (34) 

(ii) there are no two data locations totally redundant, 
i.e,, 

C(z a -Zfi) = C{x4-x p ) y fbra11lal v ..'.,n t 
if and only if : a = o/ 

Independence case 

If the n KV data Z(z a ) are independent (or uncor- 
rected) one from another, the data covariance matrix 
[C(x a — Xfi)] reduces to a diagonal matrix with all el* 
cments of the diagonal equal to the common variance 
C(0) = Var {«(*)}, for all x. 
The OK system (4) becomes: 



i.e. 



{ 



A fl C(0) + p . C(x a - x 0 ) , * « 1, . - . , n /-% 
E«A, = 1 («) 



•fijr *C(0)Il-aa-^-«o)] 

with; p(A) * C(h)/C(Q) 6 bdng the standard- 

ized covariance, or a co^^elog^ana, , ' measuring the correkr 
tion between two values Z(x) 9 Z(x + h) distant of vector 

K , moreover, the n data are uricorrelated with the un- 
known Z(x 0 ) : . 

/>(x a -a*) = 0 , or = l, — ,n, the OK system (6) then 
yields: 

A a = l/n, a = l,»-*,n 
-/i = C(0)/n 



{ 



(7) 



Thus: £•(*„) == iESi, whatever the location rr b 
and: <r 0Jr = C(0) + > 0 

In the case of total spatial independence, the OK es- 
timate reduces to the aritlmietic mean of the n data re- 
tained, 

In the SK case Lesson II (16), the SK estimate did 
reduce to the known mean mo. 1& the OK case, that 
mean is unknown and is estimated by the mean of the n 
data. 

The exactitude property ■ \ 1 ' 

If one datum, say, 2(x*0> is considered to be the un- 
known 4?(xo), then; 

C(x a < - x fi ) » C(x 0 > x tf ) , for all 0 sa 1,/ • - , n, and the 
OK system (4) is written: 

A*C(x a ~ */0 + ^C(x fl - x 0 ) + M = C(x 0 - x 0 ) , 
a = l,...,n 

The unique solution is: A^ =% A^ = 00 ^ o/,/i = 0. 
Thus, the OK estimate identifies the datum value Z(x Q <) = 
Z(x 0 ), and: 

^Sir-C(0)-A o *C(0> = 0 (8) 

As SK, the OK algorithm features the exactitude prop- 
erty: the OK surface, Z*(x 0 ) ,x 0 e Domain, honors the 
data values at data locations. 
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14 CHAPTER II. GETTING STARTED 

semivariogram models, which are then promptly converted into equivalent 
covariance models. The "variogram" programs of Chapter Ill allow the com- 
putation of covariance functions in addition to variogram functions. 

A Note on Generalized Covariances 

Generalized covariances of order k are defined as variances of differences of 
order (A + l) o the initial RF Z(u) (see, e.g. [39, 187]). The traditional 
variogram 2 7 (h), defined in relation (11.11). as the .variance of the .first-order 
difference of the RF Z(u), is associated to a generalized covariance^f order 
zero. The order zero stems from the variogram expression's filtering any zero- 
order Polynomial of the coordinates u, such as m(u)^constant, added to the 
, m °del Z(u). Similarly, a generalized covariance of order k wouloVfilter a 
polynomial trend of order fc added to the RF model Z(u). 

Unfortunately, inference of generalized covariahces of order i poses 
severe problems since experimental differences of order Jfc + 1 are not readilv 
available if the data are not gridded. In addition, more straightforward algo- 
rithms for handling polynomial trends exist, including ordinary kriging with 
moving data neighborhoods [112]. ,., 6 . 5 

tp J? 6 /^ 01 " 6 ,'. noiwi t hstandin S the possible theoretical importance iof the 
IRF-A formalism, it has been decided not to include generalized covariances 
and the related intrinsic RF models of order k in this version of QSLIB. ~ 

II. 1.4 Kriging | 

Kriging is "a collection of generalized linear regression technique!: for mini- 
mizing an estimation variance defined from a prior model for a covariance" 
([146], p. 41). ■ 

Consider the estimate of an unsampled value *(u)- from neighboring data 
values z( Ua ),a = 1 .,„. The RF model: Z(u) is stationary with mean m 
and covariance C(h). In its simplest form, also known as simple kriging (SK) 
the algorithm considers the following linear estimator! • 



Z'sk(") = £ A Q (u)Z(u Q ) + (l-£ A«(u)l m . (11.12) 



The weights A Q (u) are determined" to minimize the error variance, also called 
the estimation variance." That minimization results in a set of normal 
equations [102, 124]: . 

n 

£ x 0(u)C(up - u a ) = C(u - u a ), (H.13) 
0=1 



Va = l,...,n 
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1 tec hriK^^b r mini- 
'del fc^a cbvariance" 



"om nejg^Bb'f iri'g*data 
:ionary*Wi tB; mean rn 
^sinipl^ikri^ngtSK), 



(H12). 



variarp^Iso called 1 
in a .setiof.noniiaJ 



(11.13.) 



The corresponding minimized estimation variance, or kriging variance, is: 

n 

, vl K (u) = C(0) - £ A Q (u)C(u - n.).> 0 (n.14) 



Q = l 



, Ordinary- kngmg (OK) is the most commonly used variant of the previous 
simple kriging algorithm, whereby the sum of the weights T n _ A (u) is 
constrained to equal 1. This allows building an estimator ZA^u) 'that does 
not require prior knowledge of the stationary mean m, yet remains unbiased 
• in the sense that E {Z^ K (u)} = E {Z(u)}/ 

Non-linear kriging is but linear kriging performed on some non-linear 
. transform of the ,-data, e.g., the log-transform In z provided that z > 0 
or the indicator transform as defined in relation (II 6) 

Traditionally, kriging (SK or OK) has been performed to. provide a "best" 
linear unbiased estimate (BLUE) for unsampled values *(u), with the kriging 
-variance being used, to define Gaussian-type confidence intervals, e.g., 

• -*Prob{Z(u) e'[z' SK (u) ± 2<t sk (u)]} « 0j9S ';. 

of th p n d°m ely ' kri . ging ° f thC type (IU4 )' bein 6 independent 

confiLra IT ? ■ ° * * CaBp " iw ° f Amative geometric data 

aSatv [99} ' & ^ ***** ™* ° f ldcal estimation 

Til S^^r", ^ u°T t0 r ' aliZC that kri ging estimators of type 
(11.12) are best only m the least-squares error sense for a given covari- 
a,ce/var,ogram model. Minimizing an expected squared error need nol be 
the nipstrelevant estimation criterion for the study at hand; rather, one might 

2k?J J }" This.deasion-analysis approach to estimation requires a 
: P™baMitv d 1 str 1 but 10 noftype (II.2), Prob{Z(u) < z|(n)}, for the RV Z(u) 

„ remarks seem to imply the limited. usefulness of kriging and geo- 

tat.st.cs. as a whole. Fortunately, the kriging algorithm has two chafac- 
r, 1C prpperties that allow its use in determining posterior ccdfs of type 
•2r^^' C,6rMc Pr ° PertieS m the basis forV respectively, the 
^r^^^ - ^ «"*•* <*> aPP-h to 

(i) The Multi-Gaussian Approach: If the RF model Z(u) is multivariate 
r Russian then the simple kriging estimate (II. 12) and variance (II 14) 
identify the mean and va riance of the posterior ccdf! In addition, since 

form! «n S * mpl % hist0 « ram is " ot normal, a normal score-transform (a non-linear trans- 
iLrf ? b « P^formed on the original .-data. A multi-Gaussian model K(u) s ?h en 
hWu i 1 n ° rn ! a, 1 . 3Cbre data - For «™ P le, the transform is the logarithm L(u) - 
Mu> ,f the ,-sample histogram is approximately lognormal. Kriging and/or s7mu at on 

- a ut f s« : imed °" ^ WHh tHe reSultS appropriately bafk-transfor^ 

- values (see programs nscore and backtr in Chapter VI). 



| 'W \'}\ 1 1 || ( |i ji|r< || ii ' 



